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ABSTRACT 

We present seven light curves of the exoplanet system HAT-P-3, taken as part of a 
transit timing program using the RISE instrument on the Liverpool Telescope. The 
light curves are analysed using a Markov-Chain Monte-Carlo algorithm to update 
the parameters of the system. The inclination is found to be i = 86.75*021 degj the 
planet-star radius ratio to be R p /R+ = 0.1098* o i ' o , q 12 \ and the stellar radius to be i?* = 
0.834+q'q26^©' consistent with previous results but with a significant improvement in 
the precision. Central transit times and uncertainties for each light curve are also 
determined, and a residual permutation algorithm used as an independent check on 
the errors. The transit times are found to be consistent with a linear ephemeris, and 
a new ephemeris is calculated as T c (0) = 2454856.70118 ± 0.00018 HJD and P = 
2.899738 ± 0.000007 days. Model timing residuals are fitted to the measured timing 
residuals to place upper mass limits for a hypothetical perturbing planet as a function 
of the period ratio. These show that we have probed for planets with masses as low 
as 0.33 M® and 1.81 in the interior and exterior 2:1 resonances, respectively, 
assuming the planets are initially in circular orbits. 

Key words: methods: data analysis, stars: individual (HAT-P-3), planetary systems, 
techniques: photometric 



1 INTRODUCTION 

Transiting planets are vital to our understanding of the 
structure of planetary systems, as they allow measurement 
of the radius of the planet, and the mass of the planet un- 
ambiguously when coupled with radial velocity data. This 
allows the density and surface gravity to be derived, and the 
internal structure of the planet may be inferred. On-goin g 
transit su rveys such as Supe r WASP ((Pollacco e t al. 20081) . 
HAT Net (|Bakos et all 12002] ) and XO (|McCullough et all 
120051 ) are now pushing the number of known transiting 
planets towards 603- Most of the transiting planets are hot 
Jupiters, which produce a ~1% dip in the flux when tran- 
siting their host stars. However, an Earth-sized planet will 
produce a dip of only ~ 0.01%, and expensive space-based 
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transit surveys (e.g. C0R0T, Kepler) are required to reach 
this level of accuracy. 

Another method to detect Earth-sized planets is 
through measuring Transit Timing Variations (TTV). A 
transiting planet will maintain a constant period whilst or- 
biting its host star, excluding tidal effects and general rela- 
tivity. Therefore, if we measure the time of mid-transit for 
a given system, and see variations from a constant period, 
we may conclude that there must be a third body in the 
system perturbing the orbit of the known transiting body 
llMiralda-Escude1l2002l:lHolman fc Murray! [20051 ; I Ago! et all 
l2005l : lHevl fc Gladmanll2007l ). TTV is particularly sensitive 
when the third body is in a resonant orbit, in which case sub- 
Earth mas s planets may be detected. It is also sensitive to 
exomoons jKippindl2009l ) and Trojans (|Ford fc Gaudill2006l : 
iFord fc Holmanl I2007T ) . and therefore has the potential to 
provide the first detection of an Earth-mass body orbiting 
a main-sequence star other than our own. 
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Detecting Earth-mass bodies in resonant orbits requires 
the central transit times to be measured to an accuracy of 
~10 seconds, for which we need high precision light curves 
obtained at high cadence. In theory, we can measure tran- 
sit times to several seconds with a moderate-sized telescope, 
but in practice we are limited by correlated noise in the 
light curves, which arises from unknown changes in either 
the CCD sensitivi ty, telescope optics or observing conditions 
l|Pont et al.ll2006T ). There may also be real brightness vari- 
ations in the flux of the target or comparison stars. When 
these sources of correlated noise are minimised, it is possi- 
ble to measure central transit tim es to < 10 seconds (see 
iPont et a.1 .1120071 : IWinn et al.| [2009. for examples using space 
and ground-based photometry, respectively) . This allows the 
detection of Earth-mass planets in low-order mean-motion 
resonance, or more massive planets out of resonance (se e e.g. 
ISteffen fe Agolll2005l : IXgol fe Steffenll2007l ; lBeanll2009l '). 

We developed the RISE (Rapid Imager to Search for Ex- 
oplanets) instrument to obtain high precision light curves of 
exoplanets for TTV measurements. RISE is a fast readout 
camera mounted on the 2.0-m Liverpool Telescope (LT) on 
La Palma. It was commissioned in 2008 February, and since 
then there have been on going observations to observe exo- 
planct transits and detect a TTV signal. F i rst re sults from 



RISE have been present ed in Gibson et al. l|2008l . hereafter 



G08) for WASP-3, and in lGibson et all l|2009l , hereafter G09) 
for TrES-3, where ~10 seconds timing accuracy is achieved 

for the best transits. 

HAT-P-3 was discovered by iTorres et all l|2007l . here- 
after T07), and is a ~0.6 Mj up mass planet orbiting a K- 
type dwarf star with a period of ~ 2.9 days. The best light 
curve was taken with the 1.2 m telescope at FLWO, and 
since then there have been no published high quality light 
curves of this system. Here, we present a further seven RISE 
transit light curves of HAT-P-3, and use them to improve 
the system parameters. Additionally, we extract the TTV 
signal in an effort to detect a third body in the system. 

In Section [5] we describe the observations and data re- 
duction, and in Section [3] describe how the light curves are 
modelled and in particular how the central transit times and 
uncertainties are found. Our results are presented in Sec- 
tion^ where we use the transit timing residuals to place up- 
per mass limits on a perturbing planet that could be present 
in the HAT-P-3 system without being detected from our ob- 
servations. Finally, in Section [S] we summarise and discuss 
our results. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 RISE photometry 

Five full and two partial transits were observed using the LT 
and RISE from 2009 January 24 to 2009 May 26 . The RISE 
instrument is described in detail in ISteele et all <|2008f > and 
G08. It consists of a frame transfer CCD with a relatively 
large field-of-view (9.4 x 9.4 arcmin 2 ), and a single wide 
band filter (~ 500-700 nm). 

For all observations, the instrument was in 2 x 2 bin- 
ning mode, giving a scale of 1.1 arcsec pixel" 1 . For the first 
two transits, exposure times of 5 seconds were used with 
the telescope slightly defocused to reduce flat fielding er- 
rors and avoid saturating the CCD. A total of 2 520 images 



were obtained over a 3.5 hour period, allowing ~35-45 mins 
both before and after the transit. For the remaining tran- 
sits, exposure times of 4 seconds were used, and the tele- 
scope had slightly less defocusing, to allow 3 150 images to 
be taken over a ~3.5 hour period for the full transits. Due 
to observing constraints and weather, less than 3 150 images 
were taken on most nights, and the actual number of images 
obtained each night is shown in Table [1] The images have 
a typical FWHM of ~2-4 pixels (~2. 2-4. 4 arcsec), and the 
nights were clear for the majority of the observations, ex- 
cept for the night of 2009 April 24, where large scatter due 
to clouds can be seen in the light curve. 

Images were first de-biased and flat-fielded with com- 
bined twilight flats using standard IRAF0 routines. Aper- 
ture photometry was then performed on the target and com- 
parison stars with Pyra]f] and the DAOPHOT package. The 
number of comparison stars and aperture size for each of the 
nights is shown in Table [T] These varied as the conditions 
and field orientation changed for each night of observations, 
and were selected to minimise the out-of-transit rms. The 
light curves were then extracted by dividing the flux from 
HAT-P-3 from the sum of the flux from the comparison stars 
(all checked to be non- variable) . Initial estimates of the pho- 
tometric errors were calculated using the aperture electron 
flux, sky and read noise. The light curves were then nor- 
malised by dividing through with a linear function of time 
fitted to the out-of-transit data, setting the unocculted flux 
of HAT-P-3 equal to one. The normalisation parameters are 
allowed to vary freely when fitting the light curves, to ac- 
count for any errors resulting from this procedure. The light 
curves, along with their best fit models and residuals (see 
Section r3.1|l . are shown in Figures [1] and [2] in 1 minute bins. 
A phase folded light curve generated from the seven RISE 
light curves is shown in Figure [3j again in 1 minute bins. 



3 LIGHT CURVE MODELLING AND 
ANALYSIS 

3.1 Determination of system parameters 

In order to determine the system parameters from the tran- 
sit light curves, a parameterised model was constructed as in 
G08 and G09. This used Kepler's Laws and assumed a cir- 
cular orbit to calculate the normalised separation (z) of the 
planet and star centres as a function of time from the stel- 
lar mass and radius (M* and -R*), the planetary mass and 
radius (M p and R P ), the orbital period and inclination (P 
and i), and finally a central tr ansit time for each light curve 
(T ,„). The analytic models of iMandel fe Agoll (|2002f ) were 
then used to calculate the stellar flux occulted by the planet 
from the normalised separation and the planet/star radius 
ratio (p) assuming the quadratic limb darkening function 



a(l- M )-6(l- M ) 2 



2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

3 Pyraf is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 
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Table 1. Summary of the RISE light curves of HAT-P-3, showing the number of 
comparison stars and aperture size used, and the rms of the residuals after placing 
in 1 minute bins. 



Night 



Exposures No. comparison Aperture size rms (residuals) 











stars 


(pixels) 


(mmag) 


2009 Jan 24 


2520 


X 


5s 


5 


7 


0.88 


2009 Jan 27 


2520 


X 


5s 


5 


8 


1.27 


2009 Mar 29 


2700 


X 


4s 


2 


8 


1.37 


2009 Apr 24 


3140 


X 


4s 


2 


5 


1.99 


2009 Apr 27 


3150 


X 


4s 


3 


7 


0.96 


2009 May 23 


2851 


X 


4s 


3 


8 


1.01 


2009 May 26 


3150 


X 


4s 


2 


8 


1.44 



1.01 
1.00 
0.99 
0.98 
0.97 
1.01 
1.00 
0.99 
0.98 
0.97 
1.01 
1.00 
0.99 
0.98 
0.97 
1.01 
1.00 
0.99 
0.98 
0.97 




0.60 



0.64 



0.68 0.72 
HJD - 2454856 



0.76 




0.52 



0.56 



0.60 0.64 
HJD - 2454859 



0.68 
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HJD - 2454920 
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Figure 1. RISE light curves of HAT-P-3 taken from 2009 Jan- 
uary 24 to 2009 April 24 in 1 minute bins. Their best fit models 
are over-plotted and the residuals from the best fit are shown 
offset below each light curve. 



where / is the intensity, /j, is the cosine of the angle between 
the line-of-sight and the normal to the stellar surface, and 
a and b are the linear and quadratic limb darkening coeffi- 
cients, respectively. 

Limb darken i ng p arameters were obtained from the 
models of IClaretJ (|200fj| ). We linearly interpolated the AT- 
LAS tables for T eff = 5185K, logg = 4.61, [Fe/H] = 0.27 
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Figure 2. Same as Figure[T] for the light curves taken from 2009 
April 27 to 2009 May 26. 

and v t = 2.0 kms^ 1 (from T07) to obtain limb darkening 
parameters in both the V and R bands. The average from 
the V and R bands was then adopted as our theoretical limb 
darkening parameters. Several tests were performed to ex- 
amine the effects of the choice of limb darkening parameters 
on the results, which are described at the end of this section. 

A Markov- Chain Monte-Carlo (MCMC) algorithm was 
then used to obta in the best fit parameters and their uncer- 
tainties (see e.g., iTegmark et al J 12004 iHolman et aTl 120061 : 



ICollier Cameron et al.l 120071 : IWinn et al l |200S| ). This 
sists of calculating the \ 2 fitting statistic 



2 _ (/j.ofts — fj,calc) 

3=1 3 



where fj.ots is the flux observed at time j, aj is the corre- 
sponding uncertainty and fj yCa ic is the flux calculated from 
the model for time j and for the set of physical parameters 
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Figure 3. Phase folded light curve of the seven RISE light curves 
shown in 1 minute bins. 

described above. M*, M p and P were held fixed at the val- 
ues obtained in T07, and we fitted for i, p, R* and To, n . The 

1/3 

scaling relation R+ oc MJ was then u sed to propagate the 
error m M* to i?* l|Holman et al.ll2006t ). Subsequent param- 
eter sets are then chosen by applying small perturbations to 
the previously accepted set. New parameter sets are always 
accepted if x 2 is decreased, or with probability exp(— Ax 2 /2) 
if x 2 is increased, where A\ 2 represents the difference in y 2 
calculated for the old and new parameter sets. A further 
two free parameters are added for each light curve to ac- 
count for the normalisation, namely the out-of-transit flux 
(foot,n) and a time gradient (tc ra d,n)- The procedure is sim- 
ilar to that used in G08 and G09, to which the reader is 
referred for details. 

Photometric errors <jj were first re-scaled so that the 
best fitting model for each light curve had a reduced \ 2 of 
1. They were then re-scaled once more by a factor 8 to ac - 
count for correlated noise according to W inn et all (2008), 
again described in G08 and G09. The value determined for 
8 depends strongly on the choice of averaging times used to 
analyse the residuals. Usually the average time in a given 
range is used. However, in the this analysis we use the max- 
imum value for 8 in the range 10-35 minutes as in G09, to 
be as conservative as possible in the resulting errors. The 
values of 8 for each light curve are shown in Table [3] 

An initial MCMC analysis was used to estimate the 
starting parameters and jump functions for p, i, i?*, Tb, n , 
foot,n and tcrad,n- An MCMC run was then started for all 
seven light curves, consisting of five separate chains each 
with 500 000 points. The initial parameter set was chosen 
by adding a 5a Gaussian random to each parameters pre- 
viously determined best fit value, and the first 20% of each 
chain was eliminated to keep the initial conditions from in- 
fluencing the results. The remaining parts of the chains were 
merged to obtain the best fit values and uncertainties for 
each free parameter. The best fit value was set as the modal 
value of the probability distribution, and the la limits to the 
values where the integrals of the distribution from the mini- 
mum and maximum values were equal to 0.159. To test that 
the chains had all converged to the same region of parame- 
ter sp ace, the Gelman & Rubin statistic (jGelman fc Rubin 
1992) was then calculated for each of the free parameters, 



and was found to be less than 0.5% from unity for all pa- 
rameters, a good sign of mixing and convergence. 

To account for errors in the limb darkening parameters, 
the linear limb darkening coefficient (a) was allowed to vary 
freely whilst holding the quad ratic limb darkeni ng coefficient 
fixed at its theoretical value (|Southworthll2008l ). The chains 
converged using this technique, and did not require a prior 
to be set on the limb darkening as required for fitting the 
TrES-3 light curves in G09. This is likely due to the smaller 
impact parameter of HAT-P-3 compared to TrES-3. The 
resulting uncertainties are larger when varying the linear 
limb darkening parameter, and therefore these values are 
adopted as our final MCMC solutions. A further check on 
the limb darkening parameters was performed by adopting 
the limb darkening coefficients for the individual V and R 
filters, and repeating the above procedure, which led to no 
significant changes in our results. 



3.2 Central Transit Times 

The central transit times from the MCMC fit described 
above were adopted as the best fit values. Two different 
methods were used to account for red noise in the light 
curves, which is the biggest problem when trying to obtain 
robust central transit times. The first was re-scaling the er- 
ror bars by a factor prior to the MCMC fits. 

The second method was to use a res idual-permutation 
(RP) or "prayer b ead" algorithm (see e.g. ISouthworth||2008l ; 
iGillon et al.ll2009h . as described in G09. This consists of re- 
constructing each light curve from its best fit model and 
residuals (from the MCMC fit) by adding them, each time 
shifting the residuals by a random amount. For each recon- 
structed light curve, a fit is performed by minimising x 2 > 
and the errors are estimated from the resulting distribution 
of parameters. Twenty thousand such fits were performed 
for each light curve, with Af*, R+, i and p selected from a 
Gaussian distribution at the start of each. The values and 
uncertainties were taken from T07 for M*, and from the 
MCMC fits of the RISE light curves for i, p and i?*. T ,„, 
foot,n, Tcrad,n and a were allowed to vary freely, with their 
starting points set by randomly selecting a value within 10<r 
from the MCMC best fit values. 

For each light curve, the largest uncertainty in the cen- 
tral transit time was selected for the TTV analysis. In all 
cases this was from the MCMC fitting, which gave errors 
ranging from ~1.07-1.80 larger than the RP method. This 
is due to selecting the largest value in a given time range, 
rather than the average value. 



4 RESULTS 

4.1 System parameters 

The system parameters derived from the MCMC fits of 
the RISE light curves are given in Table [5] We find p — 
0.1098^;^, i = 86.751 ,; 2 . 2 deg and i?* = 0.834t°°^_R Q , 
consistent with the values determined in T07, but with 
smaller uncertainties. The stellar mass is assumed from 
T07, and the planet radius and density were determined as 
0.890lo °29-R./ up and 1.054+°' J^g cm" 3 , respectively. Again, 
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Table 2. Parameters and la uncertainties for HAT-P-3 as derived from 
MCMC fitting of RISE light curves and some further calculated parameters. 



Parameter 


Symbol 


Value 


Units 


Planet /star radius ratio 


P 


n inQ8+ ' 0010 
U.1U9S_ 0012 




Orbital inclination 


i 


86.75±°-j2 


deg 


Impact parameter 


b 


57R+ - 022 
u -°'°_0.033 




Transit duration 


T d 


9 n«7+ - 018 


hours 


Transit epoch 


To 


2454856.70118 ± O.OOOf 8 


HJD 


Period 


P 


2.899738 ± 0.000007 


days 


Stellar radius" 


R* 


n oq/,+0.018 

u - oo *-0 .026 


R Q 


Planet radius 


Rp 


0.890i± 


Rj 


Planet mass 6 


M v 


0.599 ±0.026 


M., 


Planet density 


Pv 


1 054+ ' 113 


g cm~ 3 


Planetary surface gravity 


log 9p 


n Q7 n + 0.()33 

J^'J_ .029 


[cgs] 



1/3 

"After error propagated according to ij* oc MJ , using M* from T07. 



b From T07, displayed here for convenience. 



Table 3. Central transit times and uncertainties for the RISE 
photometry including the error source. 



Epoch 


Central Transit Time 
[HJD] 


Uncertaintly 
(days) 


/3 C 





2454856.70137 


0.00024 


1.44 


1 


2454859.60024 


0.00037 


1.53 


22 


2454920.49567 


0.00055 


1.76 


31 


2454946.59260 


0.00065 


2.02 


32 


2454949.49334 


0.00040 


2.17 


41 


2454975.59037 


0.00034 


1.61 


42 


2454978.48993 


0.00051 


1.63 



c Re-scale factor from red noise analysis (see Section \3.1) . 



these parameters are consistent with T07, but with smaller 
uncertainties. 



6 r 

4 - 




-4 - 

_6 I > ' > ' ' >- 

-200 -150 -100 -50 50 

Epoch 

Figure 4. Timing residuals of the RISE transits (triangles) 
and the discovery epoch from T07 (square) using the updated 
ephemeris. 



4.2 Transit ephemeris 

The central transit times for the RISE light curves are shown 
in Table A new ephemeris was calculated by minimising 
X 2 through fitting a linear function of Epoch E and Period 
P to the transit times, 

T C (E) = T c (0)+EP, 

where E = was set to the transit from 2009 January 24 
taken with RISE, as it has the smallest uncertainty. We 
included the discovery epoch from T07. The results were 
T c (0) = 2454856.70118 ± 0.00018 and P = 2.899738 ± 
0.000007. Figure [3] shows a plot of the timing residuals of 
the transits using this updated ephemeris. 

For the data, a straight line fit gives \ 2 = 6.79 for 6 de- 
grees of freedom, or a reduced Xred ~ 1-13. Thus no signif- 
icant timing signal is found in the timing residuals. Indeed, 
removing all transits with less than ~20 minutes either be- 
fore ingress or after egress (Transits E = 1, 22 and 41), as 
suggested in G09, results in a \ 2 of 2.37 for 5 degrees of free- 
dom (Xred = 0.79), confirming that a straight line provides 
a very good fit to the transit times. 



4.3 Limits on a second planet in the HAT-P-3 

system. 

Whilst not showing a significant TTV signal, upper mass 
limits may still be placed on the presence of a third body in 
the system as a function of period ratio, by fitting model tim- 
ing residuals to the timing residuals measured with RISE. 
The same procedure was used in G09, where the equations 
of motion were integrated using a 4th-order Runge-Kutta 
method, with the first two bodies representing the star and 
planet of the HAT-P-3 system, and the third body repre- 
senting a hypothetical perturbing planet. The transit times 
were then extracted when the star and transiting planet 
align along the direction of observation, and the residuals 
of a straight line fit are the model timing residuals. 

Due to limitations in computation, we assume the am- 
plitude of the timing re siduals is proportional to the mass of 
the p erturbing planet l|Agol et al.l 2005; Hol man fc Murravl 
2005). We also assume that the timing residuals increase 
with increasing eccentricity of the perturbing planet. Hence, 
to set upper mass limits we can assume the planets have ini- 
tially circular orbits. This is not necessarily true near mean- 
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motion resonance (G09) , which is discussed later. The orbits 
of the planets are assumed to be coplanar throughout. 

Models were created for an Earth-mass perturbing 
planet, with a perturbing planet to transiting planet period 
ratio distributed from 0.2 to 5.0. The sampling was increased 
around both the interior and exterior 2:1 resonance (where 
we expect to probe for the smallest planets), and the transit 
times were extracted at six directions of observation starting 
along the y-axis, spaced evenly around the star. Each model 
was started with the two planets aligned along the y-axis, 
at opposing sides of the star, and we simulated 2 years of 
TTVs. The upper mass limit was determined for each model 
by scaling the mass of the perturbing planet until the y 2 of 
the model fit increased by a value A\ 2 ~ 9 (St effen Sz Agoll 
120051 ; lAgol fc Steffedl2007h from that of a linear ephemeris 
(i.e. a straight line). The \ 2 is then minimised along the 
epoch, and for a constant added to the residuals, and the 
mass of the perturbing planet is scaled again until the max- 
imum mass allowed is determined. This is repeated for each 
observation direction, and the maximum mass found is as- 
sumed as our upper mass limit for each period ratio. This 
process was then repeated twice, with the initial mass of the 
perturbing planet set to the previously determined upper 
mass limit. This was found to have little effect on the upper 
mass limits set, justifying our assumption that the timing 
residuals scale with the mass of the perturbing planet. 

Figure [5] shows a plot of the upper mass limits as a 
function of the period ratio. The solid black line is the up- 
per mass limits set from the three-body simulations, and the 
horizontal dashed line represents an Earth-mass planet. The 
results show that our data were sufficiently sensitive to probe 
for masses as small as 0.33 and 1.81 M® in the interior 
and exterior 2:1 resonance, respectively. The greater sensi- 
tivity in the interior 2:1 resonance may be explained in terms 
of the libration cycle of the timing residuals. The libration 
period of the interior resonance is ~65 orbits of HAT-P-3b, 
which allows most of the libration cycle to be covered by the 
RISE transits. However, the libration period for the exterior 
2:1 resonance is longer, ~200 orbits of HAT-P-3b, so only a 
small portion of the libration cycle is covered by the RISE 
transits. Therefore, in the exterior 2:1 resonance, lower mass 
perturbers are not as easily constrained by the transit times. 

G09 found the assumption that the smallest timing 
residuals occur when the perturbing planet has an initially 
circular orbit is not valid near mean-motion resonance. This 
was tested in the same way by creating a series of models 
around the resonances with non-zero eccentricity. A set of 
models with eccentricities ranging from to 0.23 were cre- 
ated around the exterior 2:1 resonance, and it was found 
that the amplitude of the timing residuals may drop at low 
eccentricities by as much as an order of magnitude before 
rising again at high eccentricities. This may be partly due 
to the planets getting perturbed by different amounts from 
their initial perio d ratios and overshooting the resonance 
l|Agol et al.ll2005D . resulting in different period ratios when 
only varying the initial eccentricities. A similar effect is seen 
in the interior 2:1 resonance, but with the residuals only 
dropping by about a factor of ~2. This is a similar conclu- 
sion to that in G09, and again suggests that to set true up- 
per mass limits around resonance, the effect of eccentricity 
needs to be investigated in detail. These models were then 
used to place upper mass limits as before, and we found 



that more realistic upper mass limits would be ~0.7 - 1.0 
and ~30 - 40 M® in the interior and exterior resonance, 
respectively. Some similar sets of models were produced out 
of resonance, in which the amplitude of the TTV signal does 
increase with eccentricity of the perturbing planet, and thus 
the upper mass limits found are valid. 

iBarnes fc Greenberd (|2006T ) explore the stability lim- 
its in exoplanet systems, and provide an inequality to test 
whether a system is Hill stable (equation 2). Using this 
inequality for the HAT-P-3 system (assuming an Earth- 
massed perturber), places lower and upper limits on the pe- 
riod ratio of 0.74 and 1.37, respectively. The resulting region 
not guaranteed to be Hill stable is marked on Figure [5] by 
the grey shading, although stable configurations may still 
occur in this region. Trojan compa nions may also exist in 
stable orbits near the 1:1 resonance. iMadhusudhan &; Winnl 
(2009) placed a 2a upper mass limit of ~ 260 Earth masses 
on a Trojan companion to HAT-P-3b by combining transit 
observations and radial velocity data. 



5 SUMMARY AND DISCUSSION 

This paper presents the second Transit Timing study using 
the RISE instrument on the LT, consisting of seven light 
curves of the exoplanet system HAT-P-3. An MCMC al- 
gorithm was used to determine the light curve parameters, 
which were found to be in agreement with the discovery pa- 
per, but with a higher precision. 

The central transit times and uncertainties were also 
determined using the MCMC algorithm, and a residual per- 
mutation algorithm was used as an independent check on 
the errors. Uncertainties in the transit times range from 
20 — 60 seconds, much greater than the typical transit times 
obtained for TrES-3 using RISE (G09). This is partly due to 
the shallower transit of HAT-P-3, but also due to the higher 
levels of systematic noise in the light curves. 

The central transit times were found to be consistent 
with a linear ephemeris, with Xred < 1 after removing those 
light curves with less than 20 minutes of data either before 
ingress or after egress. These transit times were then used 
to place upper mass limits of an additional planet in the 
system that could perturb the orbit of HAT-P-3b yet not be 
detected by our measurements. This showed that we probed 
for masses as low as 0.33 M® and 1.81 M® in the interior and 
exterior 2:1 resonance, respectively. However, larger planets 
may exist in low eccentricity orbits in the 2:1 resonances, 
and further investigation is required to explore true upper 
mass limits in a higher dimensional parameter space. These 
are comparable with the mass limits set using the same tech- 
niques for TrES-3 in G09, which had many more light curves 
available, plus transit times with smaller uncertainties. This 
is because the period of HAT-P-3 is larger than that of TrES- 
3, and reflects that the sensitivity of the Transit Timing 
method increases with the period of the transiting planet. 
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Figure 5. Upper mass limits of a hypothetical 2nd planet in the HAT-P-3 system as a function of period 
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